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Abstract 


Coronaviruses express a de-ubiquitinating protein, the papain-like protease-2 (PLP2), that re¬ 
moves both ubiquitin and the ubiquitin-like Interferon (IFN) Stimulated Gene 15 (ISG15) protein 
from target proteins. ISG15 has antiviral activity against a number of viruses therefore, we exam¬ 
ined the effect of ISG15 conjugation (ISGylation) in a model of acute viral hepatitis induced by 
the murine hepatitis virus (MHV)-3 coronavirus. Mice deficient in the ISG15 deconjugating en¬ 
zyme, ubiquitin specific peptidase-18 (USP18), accumulate high levels of ISG15-conjugated pro¬ 
teins and are hypersensitive to type I IFN. Infecting US P18"'" mice with MFFV-3 resulted in ex¬ 
tended survival (8 ± 1.2 vs. 4 days), and improved liver histology, a decreased inflammatory re¬ 
sponse, and 1-2 logs lower viral titers compared to USP18 +/ mice. The suppression of viral rep¬ 
lication was not due to increased IFN, since infected USP18" " mice had neither increased hepatic 
IFN-a, -(3 or -y mRNA nor circulating protein. Instead, delayed MFFV-3 replication coincided 
with high levels of cellular ISGylation. Decreasing ISGylation by knockdown of the ISG15 El 
enzyme, UbelL, in primary USP18 /+ and USP18" " hepatocytes led to increased MHV-3 replica¬ 
tion. Both in vitro and in vivo, increasing MHV-3 titers were coincident with increased PLP2 
mRNA and decreased ISGylation over the course of infection. The pharmacologic inhibition of 
the PLP2 enzyme in vitro led to decreased MHV-3 replication. Overall, these results demonstrate 
the antiviral effect of ISGylation in an in vivo model of coronavirus-induced mouse hepatitis and 
illustrate that PLP2 manipulates the host innate immune response through the ISG15/USP18 
pathway. 
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45 Statement of Importance 

46 There have been a number of serious worldwide pandemics due to widespread infections by 

47 Coronavirus. This virus (in its many forms) is difficult to treat, in part because it is very good at 

48 finding "holes" in the way that the host (the infected individual) tries to control and eliminate the 

49 virus. In this study we demonstrate that an important host viral defence - the ISG15 pathway - is 

50 only partially effective in controlling severe Coronavirus infection. Activation of the pathway is 

51 very good at suppressing viral production, but over time the virus overwhelms the host response 

52 and the effects of the ISG15 pathway. This data provides insight into the host-viral interactions 

53 during Coronavirus infection and suggests that the ISG15 pathway is a reasonable target for con- 

54 trolling severe Coronavirus infection, though the best treatment will likely involve multiple 

55 pathways and targets. 
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Introduction 


Coronaviruses cause both common and severe clinical illness, as manifested by the Se¬ 
vere Acute Respiratory Syndrome (SARS) epidemic. During the 2002-2003 SARS epidemic, 
8422 people were affected of whom 916 died from acute respiratory distress syndrome (1,2). 

The episodic re-emergence of severe coronavirus infections, most recently in the fall of 2012 and 
continuing into 2013 highlights the need for treatments against Coronavirus infections (3). 

Coronaviruses are positive stranded enveloped RNA viruses with some of the largest vi¬ 
ral genomes ranging between 26-32kDa. Coronaviruses target a myriad of distinct animal hosts 
and cause disease in a number of organs such as the brain, liver and lung. Organ damage from 
severe coronavirus infections is generally the result of an over-exuberant activation of host in¬ 
nate immune mechanisms (4, 5). For example, Murine Hepatitis strain-3 (MF1V-3) infection of 
susceptible mouse strains is a model of strong innate immune activation, resulting in fulminant 
viral hepatitis (6). Following infection by MHV-3, mice die in 3-4 days of hepatic parenchymal 
destruction mediated by a robust activation of local innate immunity (7-9). Understanding the 
limits of host immunity in the MF1V-3 model may identify novel targets for the treatment of se¬ 
vere Coronavirus infections. 

Interferon (IFN) stimulation as well as bacterial and viral infection induces the expression 
of IFN stimulated genes (ISGs). One of the most abundantly expressed ISGs is ISG15, a 15-kDa 
protein (10). ISG15 is conjugated to target proteins, the process of ISGylation, through consecu¬ 
tive interactions with an El activating enzyme (UbelL), an E2 conjugating enzyme (UbcF16 or 
UbcF18) and an E3 ligase (Flerc5, EFP, FIFIARI) (11). Over 160 proteins have been identified as 
targets of ISG15 including proteins involved in processes such as protein translation, cell cycle 
regulation and immune regulation (11-13). Several ISG15 downstream targets are involved in the 
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regulation of IFN signalling, including retinoic acid-inducible gene 1 (RIG-1), signal transducer 
and activator of transcription -1 (STAT-1), janus activated kinase-1 (JAK-1) and MxA (12-14). 
Although the ISG15 protease function has the potential to modulate IFN responses, this role has 
not been consistently demonstrated (15). The effect of ISGylation on downstream proteins is 
multifaceted as ISGylation has been reported to disrupt target protein function and/or alter cellu¬ 
lar localization (16, 17). 

The conjugation of ISG15 to protein targets is offset by the deconjugating activity of the 
ubiquitin specific peptidase-18 (USP18). USP18, like ISG15, is an ISG, which is up-regulated 
after stimulation by IFN, or by viral and bacterial infection. USP18 is specific for ISG 15 and 
strips ISG 15 from target proteins through its isopeptidase activity (18). USP18" mice have 
markedly increased cellular ISGylation, and are hyper-sensitive to the effects of Type 1 IFN (19- 
21). The latter effect is independent of the isopeptidase activity, and instead due to the ability of 
USP 18 to modulate IFN signalling by binding to the IFNa receptor 2 (IFNAR2) (22). Thus, 

USP 18 has several functions important to the host innate immune response: by binding to the 
IFNAR2 it can modulate the IFN response; and, through its ISG 15 isopeptidase function it regu¬ 
lates the cellular ISGylation levels. However the implications of the balance of ISGylation and 
USP 18 isopeptidase activity require further elucidation. 

The role of ISGylation in viral lifecycles is specific to the virus involved. ISGylation can 
exert an antiviral pressure against some infections, but can also stimulate viral replication in oth¬ 
ers. For example, USP18 deficient mice have increased ISGylation and are more resistant to 
lymphocytic choriomeningitis virus (LCMV) and herpes simplex virus (HSV) models of fatal 
viral encephalitis (20). During human immunodeficiency virus (HIV) infection the conjugation 
of ISG 15 to the HIV Gag protein arrests assembly of the Gag particle on the plasma membrane 
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(23), and inhibits viral replication. On the other hand, we have found that ISGylation is necessary 
for robust production of HCV in human hepatocytes (24), and both ISG15 and USP18 are up- 
regulated in the hepatocytes of patients with chronic HCV who do not respond to exogenous IFN 
treatment (25-27). 

Previously, we demonstrated that coronavirus replication in vivo is held in check by host 
ubiquitination (28). Inhibition of the cellular proteasome leads to increased cellular ubiquitina- 
tion levels and an early interruption in coronavirus replication. Coronaviruses have evolved 
counter-measures to such host cellular antiviral mechanisms. In this case, the de-ubiquitination 
protein papain-like protease-2 (PLP2) strips ubiquitin from target proteins (29). The PLP2 pro¬ 
tein is not specific to ubiquitin: it acts on both ubiquitin and the ubiquitin-like protein, ISG15. 
This suggests that ISG15 and its conjugation to cellular proteins may also exert an antiviral pres¬ 
sure effect against coronavirus infection. Moreover, several reports indicate that PLP2 has a 60- 
fold higher affinity for ISG15-conjugated rather than Ub-modified substrates (18, 30). However, 
it is unknown whether targeting PLP2 activity in an in vitro cell culture system or in vivo will 
inhibit coronavirus replication. Furthermore, questions remain as to whether ISGylation itself is 
antiviral against coronavirus infection and/or whether the PLP2 machinery in vitro targets ISGy¬ 
lation directly to evade the host response. 

In the present study, we infected USP18" ’ mice with the coronavirus MHV-3 to evaluate 
the involvement of the ISG15/USP18 pathway to the virulence of MHV-3-induced hepatitis. We 
found that USP18" " mice are more resistant to MHV-3 infection, but that the virus gradually 
overcomes this protective effect. IFN type 1 and type 2 expression levels were not increased in 
USP18" " mice following infection by MHV-3, allowing us to study the role of increased baseline 
ISGylation in the absence of IFN signalling. Silencing ISGylation reverses the antiviral milieu 
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126 and leads to increased MHV-3 replication. Both in vitro and in vivo, viral persistence is accom- 

127 panied by increased expression of the PLP2 protein. PLP2 expression is important for allowing 

128 viral replication; specific PLP inhibitors decreased viral protein expression. Overall, these results 

129 demonstrate both the role and the limits of the antiviral effect of ISGylation in severe coronavi- 

130 rus infection. 
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Materials and Methods: 


Animals 

C57BL/6 USP18 and USP18 _/ "mice between 6 -8 weeks old were used for experiments. These 
mice were the kind gift of Dong Er Zhang (Scripps/UCSD). Animals were housed in SPF condi¬ 
tions at the MaRS-TMDT Animal Resource Centre (Toronto) and were given chow and water ad 
libitum. Following MF1V-3 infection animals were housed in sterile cages in a level 2 facility in 
the MBRC (Toronto). Mice infected with MF1V-3 were sacrificed at humane endpoints. Animals 
were treated according to the guidelines of the Canadian Council on Animal Care and were ap¬ 
proved by the University Flealth Network Animal Care Committee. 

Isolation of Peritoneal exudative macrophages (PEM) and Hepatocytes 

The isolation of PEM has been previously described (28). Briefly, mice were injected with 2ml 
of sterile 3% thioglycollate. Animals were sacrificed after 3 days and PEM were retrieved by 
washing the intraperitoneal cavity with 10ml ice cold Flank's buffered salt solution (FIBSS) (Life 
Technologies). Cells were washed 2-fold, spun down at 300 * g and re-suspended in Dulbecco's 
Modified Eagle Medium (DMEM) (Life Technologies ) supplemented with L-Gln. Using this 
method, the purity of PEM was found to be > 90% with a viability of >97% by Trypan blue ex¬ 
clusion. For experiments, PEM were plated on polystyrene plates at a density of 1 x 10 6 cells/ml 
and incubated for 24h at 37°C and 5% CCb. Non-adherent cells were then washed away. 

Primary hepatocytes were isolated as previously described (31, 32). Briefly, mice were anesthe- 
tised with 50 mg/kg of pentobarbital i.p.. The portal vein was canulated with a 21-gauge needle. 
The liver was flushed via the portal vena cava with perfusion solution (2.5mM EGTA in FIBSS 
without Ca 2+ or Mg 2+ ) at 37°C and a rate of 7ml/min using an infusion pump for 3min. The liver 
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was then perfused with solution #2 (0.02% collagenase IV (Sigma-Aldrich) with calcium and 
magnesium in HBSS) at the same pressure. The liver was then carefully harvested into a Petri 
dish containing DMEM-15 and 1 x 10" 7 M insulin (Sigma-Aldrich) and minced. This solution 
was filtered through a 120 pm nylon mesh and centrifuged for 2min at 40 x g at room tempera¬ 
ture (RT). Cells were washed 2 times with DMEM + insulin. Cells were counted on a hemocy- 
tometer and viability was determined by Trypan blue exclusion. Hepatocytes were plated at a 
density of 1 x 10 5 cells/ml in DMEM supplemented with 1 x 10' 7 M insulin on a 6-well plate 
(Corning Inc.). Hepatocytes were incubated for 2 hours and the medium was replaced with serum 
free DMEM supplemented with 1 x 10" 7 M insulin. Hepatocytes were inoculated with MHV-3 
(multiplicity of infection (MOI) = 1, 0.1 or 0.01) and incubated for lh. Cells were then washed 
with serum free medium + 1 x 10’ 7 M insulin and cultured for the indicated time. 

Virus, viral infection and viral titering 

MHV-3 virus was grown to titers of 10-50 x 10 6 PFU/ml RPMI on confluent 17CL cells. To de¬ 
termine viral titers, harvested cells or liver tissue, were lysed or homogenized in ice-cold 
DMEM-10 using a TissueLyser (QIAGEN). MHV-3 titers were assessed on monolayers of L2 
cells in a standard plaque assay (6). For in vivo studies mice were infected i.p. with 50 PFU of 
MHV-3 and housed in sterile conditions. Animals were monitored daily and sacrificed at a hu¬ 
mane end point. 

Histology 

Tissues were harvested and preserved in 10% buffered formalin. Fixed tissues were paraffin em¬ 
bedded and cut into 5pm sections by the Department of Pathology at the Hospital for Sick Chil- 
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dren (Toronto). Sections were stained with Hematoxylin and Eosin using standard protocols (33). 
Sections were scored in a blinded fashion by an experienced liver pathologist (M.J.P.). Sections 
were assessed for inflammation, parenchymal changes and tissue necrosis. 

Measurement of alanine transaminase (ALT) and (aspartate aminotransferase)AST in the 
serum 

Blood was harvested by cardiac puncture and was incubated for lOmin at RT. Samples were 
spun down at 2,000 x g and the serum was collected for analysis. ALT/AST levels were meas¬ 
ured using the Vitro DT60 II Chemistry System (Ortho Clinical Diagnostics, Neckargemund, 
Germany). 

Serum protein levels of type I and type IIIFN 

Blood was harvested from mice by cardiac puncture. Serum was incubated at 4°C overnight 
(O/N). The samples were centrifuged at 2,000 x g for lOmin and the serum was removed to a 
fresh tube. ELISA kits for IFN-a and IFN-P were purchased from PBL Interferon Source (Pis- 
cataway Township, NJ, USA). The ELISA kit for IFN-y was purchased from BioLegend Inc 
(San Diego, CA, USA). ELISA tests were performed according to manufacturer’s instructions. 

Real Time PCR 

Mice were infected with MHV-3 and sacrificed on the indicated day(s). Total RNA was ex¬ 
tracted from liver tissue using TRIzol reagent (Invitrogen) according to manufacturer’s specifica¬ 
tions. The isolated RNA was treated with DNase (Qiagen) and the quality of the RNA was as¬ 
sessed by measuring 260/280 on a Nanodrop spectrophotometer, lpg of RNA was reverse tran- 
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scribed using the First Strand cDNA Synthesis Kit (GE Flealth Care) using random hexamer oli¬ 
gonucleotides. The reverse transcriptase reaction was run as per manufacturer’s specifications. 
Primers used for Real Time PCR reactions were based on published and primer bank sequences. 
The following primer sequences were used: mouse IFN-a forward primer 5'- 
CTTTGGATTCCCGCAGGA-3'; mouse IFN-a reverse primer 5'- TGTAG- 
GACAGGGATGGCTTGA-3'; mouse IFN-P forward primer 5'- 
TGAATGGAAAGATCAACCTCACCTA-3mouse IFN-P reverse primer 5'- 
CTCTTCTGCATCTTCTCCGTCA-3'; mouse IFN-y forward primer 5'- 
CAGCAACAGCAAGGCGAAA-3'; mouse IFN-y reverse primer 5'- 

CTGGACCTGTGGGTTGTTGAC-3'; mouse hypoxanthine phosphoribosyltransferase (HPRT) 
forward primer 5'-GTTGGATACAGGCCAGACTTTGTT G-3'; mouse HPRT reverse primer 5'- 
GATTCAACTTGCGCTCATCTTAGGC-3'; mouse glyceraldehyde 3-phosphate dehydrogenase 
(GAPDH) forward primer 5'-ACAACTTTGGCATTGTGGAA-3'; mouse GAPDH reverse 
primer 5'-GATGCAGGGATGATGTTCTG-3'; MHV-3 PLP2 forward primer 5'- 
AAATGTGGCTTGTTTGATGC-3and, MHV-3 PLP2 reverse primer 5'- 
CCTT CGCT AAGACC AAGG AC-3'. 

Real Time PCR was performed on an ABI PRISM 7900H (Applied Biosystems) with SYBR 
Green real time PCR master mix (Applied Biosystems) according to manufacturer’s specifica¬ 
tions using the 96-well plate format. Data was analyzed using SDS 2.2 Software (Applied Bio¬ 
systems) using the standard curve method. Samples were normalized to HPRT and GAPDH. 

Western Blot (WB) 
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Harvested tissues were homogenized and lysed in ice-cold cell lysis buffer. Cells from in vitro 
cultures were lysed with 2 x Laemmli buffer and O.lmM DTT. Protein concentrations were de¬ 
termined by the Bradford assay. 20 gg of proteins was visualized by 10% SDS-PAGE gel and 
transferred to a nitrocellulose membrane (Pall). Membranes were probed with the following an¬ 
tibodies: anti-nucleocapsid (N) protein mouse monoclonal (made in-house), anti-mouse P-actin 
mouse monoclonal antibody (Sigma), anti-mouse ISG15 rabbit polyclonal antibody, anti-mouse 
UbelL goat polyclonal IgG, or anti-mouse STAT1 (p84/p91) rabbit polyclonal (Santa Cruz Bio¬ 
technology). The corresponding horseradish peroxidase (HRP) conjugated secondary antibodies 
were used: Sheep anti-mouse IgG (GE Healthcare), donkey anti-rabbit IgG (Santa Cruz Biotech¬ 
nology) or a donkey anti goat IgG (Santa Cruz Biotechnology). Blots were developed using the 
ECL system (Amersham Pharmacia Biotech). 

Silencing UbelL expression 

Isolated primary hepatocytes were plated at a density of lx 10 5 cells/ml in DMEM. Cells were 
incubated O/N at 37°C, 0.5% CO 2 whereupon cells reached 60-80% confluency. Cells were 
transfected with UbelL siRNA (Santa Cruz Biotechnology) as per manufacturer’s instructions. 
Briefly, cells were washed with transfection medium and then overlaid with transfection reagent 
solution and incubated for 6 h at 37°C. DMEM medium including 10% FBS was added and the 
cells were incubated a total of 24h. The medium was replaced with normal growth medium and 
cells were infected with MHV-3 (MOI=l) ± IFN-a (PBL Interferon Source). Samples were har¬ 
vested 12h p.i. using 200pl protein lysis buffer (BD bioscience), run on a 10% SDS PAGE gel 
and visualized by WB. 
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PLP2 Inhibition 


PLP2 inhibitors were generated in house according to Ratia et al (34). Isolated PEM were incu¬ 
bated in the presence of MHV-3 (MOI=l) for lh. Medium was washed away and fresh medium 
± lOOpM inhibitor 6, 100pM inhibitor 7, 100pM GRL0617 or lOOIU/ml IFN-a was added (PBL 
Interferon Source). Cells were harvested at 9h p.i., run on a 10% SDS PAGE gel and visualized 
by WB. 

Statistics 

Statistical analysis was performed using PRISM software version 5 (GraphPad Software Inc). 
Survival curves were calculated using the Log-rank (Mantel-Cox) test. 

Data analysis was performed using a 2-way ANOVA followed by a Bonferroni post-hoc test 
unless otherwise noted. Data is represented as either mean ± SD or mean ± SEM as indicated. 
Results with a p value below 0.05 were considered significant. 
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Results 


USP18' ' mice have increased survival following MHV-3 infection 

To test our hypothesis that ISGylation has an antiviral effect on murine hepatitis virus 
(MHV)-3 infection we infected USP18" " mice, which have high baseline levels of ISGylation, 
and control USP18 +/+ mice with 50 pfu. MHV-3 (20). The infection of susceptible C57BL/6 mice 
with MHV-3 induces a fulminant hepatitis with hepatocellular necrosis, sinusoid thrombosis, 
strong cytokine induction and leads to death within 3-4 days p.i. (6). USP18 +/+ (C57BL/6) mice 
survived a median of 4 days and observed the rapid development of symptoms characteristic of 
MHV-3 infection (Figure 1A). By contrast, USP18" " mice (C57BL/6) survived up to 10 days p.i., 
with a median survival of 8 days. The lengthened survival was accompanied by improved liver 
histology at all time points. On day 3, all USP18 + ' liver histology showed 40 to 50% focal ne¬ 
crosis, and by day 4, 70-90% of the liver showed evidence of confluent necrosis (Figure IB). In 
comparison, USP18" " livers on day 4 showed only small areas of focal necrosis affecting 5-10% 
of the liver. MHV-3 infection of USP18" " mice did gradually lead to widespread liver damage 
(40-50% focal necrosis), though never as pronounced as in USP18 +/ mice. Improved liver his¬ 
tology in USP18" " mice coincided with a decrease in the levels of serum transaminases (markers 
of liver injury) and decreased viral titers compared to USP18 +/+ mice (Figure 1C-E). Histology, 
viral titers and liver function in USP18" " was improved at all time points compared with 
USP18 +/+ mice. However, at day 7 p.i., the USP18 _/ " pathogenesis approached levels comparable 
to USP18 /+ mice at day 3 p.i. in all 3 of the measured criteria including histology, serum 
transaminase levels and viral titers (Figure IE). Therefore, the absence of USP18 confers an 
early antiviral advantage against MHV-3 infection. However, USP18’ " mice ultimately suc¬ 
cumbed to infection. 
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MHV-3 infection does not induce type IIFN in USP18' ' mice 

USP18" " mice have increased baseline cellular ISGylation levels and they have increased 
sensitivity to Type I IFN signalling (20). Either of the two aspects of the USP18" " phenotype 
may explain the suppression of viral replication and decreased inflammatory response seen in 
response to MHV-3 infection. In order to determine whether increased sensitivity to type 1 IFN 
contributed to the suppression of MHV-3 replication in USP18" " mice, we measured Type I IFN 
(IFNa, IFN|3) induction by MHV-3 infection. In parallel we measured Type II IFN (IFNy), as a 
control pro-inflammatory cytokine. MHV-3 infection of USP18 /+ mice induced strong mRNA 
and protein expression of IFNa, IFNp, and IFNy in the liver and serum respectively. In USP18" /_ 
mice, there was no induction of IFNa, IFNp, or IFNy hepatic mRNA and protein levels in re¬ 
sponse to MHV-3 infection at all time points studied (Figure 2A-D). Similarly, MHV-3 infection 
did not induce type I or type II IFN induction in USP18" " primary hepatocytes (data not shown). 
Finally, USP18" " peritoneal exudative macrophages (PEM) that were infected with MHV-3 did 
not show increased Type I IFN signalling, as seen by the absence of STAT1 phosphorylation 
(data not shown). These data imply that the increased antiviral activity against MHV-3 observed 
in the USP18’ ’ mice is not due to increased IFN production or signalling. 

Cellular ISGylation inhibits MHV-3 Replication 

Since the decrease in MHV-3 replication and mortality following infection of USP18" " 
mice is not due to increased IFN or Type I IFN signalling, we asked whether decreased MHV-3 
replication and mortality was due to increased levels of ISGylation (21). First, we examined 
whether the decreased proliferation of MHV-3 in USP18 in vivo could be replicated in vitro in 
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primary hepatocytes. Hepatocytes isolated from USP18 +/+ and USP18’ " mice were infected with 
MHV-3 (multiplicity of infection (MOI)=l) and viral titers were measured up to 48h p.i. (Figure 
3A). USP18 +/ " hepatocytes reached peak MHV-3 titers, 2.09 x 10 7 ± 1.29 x 10 6 PFU, between 
12-18h p.i. Viral titers remained elevated until 36h p.i, after which the levels started to decrease, 
presumably due to cell death. Viral titers increased in USP18" " hepatocytes with slower kinetics 
than USP18 +/+ mice, with nearly a log decrease in MHV-3 titers at almost all time points. How¬ 
ever, USP18" 7 ' mice reached peak viral titers at 36h at 1.03 x 10 7 ±2.31 x 10 5 PFU, more than 
18h after USP18 +/+ peak titers. These in vitro results support the in vivo finding that MHV-3 viral 
titers are suppressed in the presence of increased ISGylation. Although the viral titers in USP18" " 
hepatocytes are suppressed initially, the virus eventually overwhelms the antiviral response. This 
is similar to the extended survival seen in USP18" " mice, but these mice ultimately succumb to 
MHV-3 infection. Of note, at early time points following MHV-3 inoculation, namely lh, 3h and 
6h, viral titers of USP18 +/+ and USP18 _/ " mice do not differ significantly (2-way ANOVA with 
Bonferroni post-hoc test). This suggests that viral binding and entry into the cell may not be af¬ 
fected by interferon stimulated gene (ISG)-15 antiviral activity however, further studies are re¬ 
quired to differentiate the effects of ISGylation on viral binding and entry and viral replication. 

The amount of ISG15-conjugated proteins is markedly higher in USP18’ " hepatocytes 
than USP18 +/+ hepatocytes at all time points following MHV-3 infection (Figure 3B). The ex¬ 
pression of viral nucleocapsid (N) protein is detectable by WB starting at 6h p.i. in USP18 +/+ 
hepatocytes and continues to accumulate until 18h (Figure 3B) when peakN protein expression 
coincides with peak viral titers (Figure 3A). By contrast, in USP18" " hepatocytes, overall N pro¬ 
tein accumulation is decreased and is delayed compared to USP18 +/+ hepatocytes. These in vitro 
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data mirror the in vivo results, and suggest that high ISGylation levels are an intrinsic antiviral 
force against coronavirus in USP18" ’ mice. 

To determine whether ISGylation is causally linked to MHV-3 virus production we 
knocked down the ISG15 pathway El enzyme UbelL to block ISGylation. UbelL is the unique 
activating enzyme of the ISG15/USP18 pathway, and forms an ATP-dependent thioester bond 
with ISG15 (11, 35). UbelL knockdown in primary USP18 +/+ murine hepatocytes decreased 
UbelL protein and attenuated the increase in ISGylation that follows exposure of cells to IFN-a 
(Figure 4). UbelL silencing in USP18 + + primary cells infected with MHV-3 abrogates ISGylation 
and leads to increased N protein expression. As shown previously (Figure 3B), N protein levels 
are markedly decreased in USP18‘ /_ compared to USP18 cells at 12h p.i.. The effect of UbelL 
knockdown is more prominent in USP18’ ’ primary cells UbelL silencing markedly decreased 
protein ISGylation in USP18" " primary cells from high baseline levels of ISGylation, although 
ISGylation levels remained at elevated levels following MHV-3 infection compared to control 
hepatocytes. Following both MHV-3 infection and UbelL knock down, the levels of ISGylated 
proteins drops further compared to UbelL knock down alone (Figure 4) and, N protein expres¬ 
sion, the surrogate of MHV-3 replication, reached levels of expression higher than induced by 
MHV-3 infection alone. MHV-3 infection decreased ISGylation in both USP18 +/+ and USP18 _/ " 
mice. The former however was only observed following treatment with IFN-a which strongly 
induces ISGylation. In the absence of IFN-a, no ISGylation is detected by WB. These data sug¬ 
gest that the UbelL conjugation of ISG15 to target proteins, and the subsequent accumulation of 
ISGylated proteins mediates an antiviral effect against MHV-3. 
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The effect of PLP inhibitors on MHV-3 infection 


Since there is an observed decrease in ISGylation in USP18 _/ " following MHV-3 infection, we 
next asked whether the MHV-3 PLP2 protein is the mechanism through which MHV-3 ulti¬ 
mately escapes the effect of host ISGylation. To this end, we treated MHV-3 infected cells with 
synthetic PLP inhibitors. The PLP inhibitors are non-covalent cysteine inhibitor compounds gen¬ 
erated against the catalytic residues of the SARS Coronavirus PLpro (34). PLP inhibitors in¬ 
creased ISGylation levels in USP18 +/+ PEM. No observable change in ISGylation was observed 
in USP18" " PEM even in less saturated blots (data not shown) due perhaps to already high base¬ 
line levels of ISGylation in USP18" " PEM. USP18 +/+ PEM infected with MHV-3 showed strong 
expression of N protein at 9h p.i. Treatment of USP18 +/+ PEM with PLpro inhibitors resulted in 
decreased levels ofN protein production (Figure 5). Treatment with GRL0617, the PLP inhibitor 
with the highest binding capacity (34), resulted in the greatest decrease in viral N protein produc¬ 
tion. In MHV-3 infected USP18PEM, where baseline levels of ISGylated proteins are high 
(21), the expression of N protein is significantly decreased compared to wildtype PEM. The ex¬ 
pression of N protein is even more inhibited following treatment with all PLpro inhibitors. By 
WB analysis, GRL0617 completely inhibited the expression of N protein. Consistent with these 
data, viral titers were decreased in the presence of PLP2 inhibitors in both USP18 +,+ and USP18’ 

" PEM (data not shown). Similar data was obtained for hepatocytes (data not shown). Overall 
these results suggest that MHV-3 proliferation and evasion from the cellular antiviral ISGylation 
milieu is dependent on PLP2 activity. 

ISGylation delays onset of MHV-3 production and death in vivo 
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377 Above, we showed that ISGylation in vitro exerts an antiviral pressure on MHV-3 coro- 

378 navirus replication. However, this antiviral pressure is overcome as the virus replicates and ex- 

379 presses the deubiquitinase protein (DUB) PLP2. In order to test whether this holds true in vivo, 

380 we inoculated USP18 +/+ and USP18 _/ " animals with 50 PFU of MHV-3 i.p. and measured ISGyla- 

381 tion and PLP2 levels in whole liver tissue. As before, USP18 /+ mice became moribund and were 

382 euthanized by day 4 p.i. Liver protein ISGylation increased and peaked at day 3p.i. (Figure 6); N 

383 protein was first detectable at day 2 and increased steadily until the animals were euthanized on 

384 day 4. By contrast, USP18 _/ " mice have higher baseline liver ISGylation levels throughout the 

385 entire course of infection, and MHV-3 replication as seen by N protein expression was first de- 

386 tected at day 6 p.i.. Concurrently, we measured viral expression of plp2 mRNA. Pip2 mRNA ex- 

387 pression increased as N protein expression increased and plp2 expression was delayed in USP18" 

388 ' mice compared to control mice (Figure 7). These data are consistent with the ability of in- 

389 creased PLP2 to gradually overwhelm the antiviral effect of ISGylation in vivo by stripping 

390 ISG15 from target proteins. 
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Discussion 


During the innate immune response, ISGylation is an important barrier against viral in¬ 
fection. ISG are up-regulated following the activation of type IIFN signalling and have various 
functions in the innate immune response including the feed-back regulation of the host immune 
response and antiviral activities (36-41). ISG15 is one of the most abundantly up-regulated ISG. 
In this study we show that ISG 15 conjugation (ISGylation) has antiviral activity in a model of 
severe coronavirus infection and that this antiviral activity is independent of type I IFN signal¬ 
ling. 

Studies have shown that ISGylation is antiviral against several viruses including, influ¬ 
enza B virus, sindbis virus, sendai virus and vaccinia virus (15, 41-44). ISG 15 and IFNaP recep¬ 
tor double deficient mice, which lack ISG 15 and have an impaired ability to induce ISG follow¬ 
ing IFNa/p stimulation, were rescued against Sindbis virus induced lethality by expressing 
ISG15 (42, 43). The antiviral activity associated with ISG15 is dependent on the conjugation of 
ISG 15 to target proteins, as the mutation of the ISG 15 C-terminal isopeptidase target motif ren¬ 
dered ISGylation ineffective against Sindbis Virus (42). Flowever, ISGylation is permissive for 
certain viruses. We and others have demonstrated that cellular ISGylation is necessary for the 
efficient replication of Flepatitis C Virus (25, 45). Our data demonstrate that ISGylation is antivi¬ 
ral to MF1V-3 infection. In USP18" " mice, where ISGylation levels are high, there was a signifi¬ 
cant decrease in viral N protein expression in vitro and in vivo that coincided with prolonged 
survival. We attribute the prolongation in survival to the antiviral effect of ISGylation, since we 
did not observe the induction of a type I or type II IFN response in the USP18" " animals. The 
role of ISGylation was confirmed by UbelL knockdown experiments, in which silencing UbelL 
in primary hepatocytes reversed the antiviral activity of ISGylation. (Figure 4). These data are in 
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accordance with studies in which both UbelL" " mice, which lack the ability to conjugate ISG15 
to target proteins, and ISG15’ " mice were more sensitive to Sindbis Vims infection (43, 46) than 
wildtype animals. Therefore, our data provides another example of the antiviral nature of ISGy- 
lation activity but in a novel context, coronavirus infection. 

Coronaviruses evade the host innate immune response by interfering with the induction 
of IFN (47-52). MHV-3 infection of USP18 +/+ mice resulted in markedly increased Type IIFN 
mRNA and protein (Fig 2A-D). By contrast, MF1V-3 infection of USP18’ ’ mice resulted in little, 
if any, IFN production at either the mRNA or protein levels. At first glance this is surprising, 
since USP18" ’ mice are hypersensitive to IFN, and we expected that any IFN-(i produced upon 
exposure to MHV-3 virus would stimulate IFN-a production in a positive feedback loop (53). 
Furthermore, the increased ISGylation observed in USP18’ " cells would be expected to lead to 
increased IFN production since the direct ISGylation of interferon regulatory transcription fac¬ 
tor-3 (IRF-3) prolongs the activation of IRF-3 and results in increased expression of IFN-p and 
other ISGs (41). However, we observed an almost complete absence of type I IFN production in 
response to MHV-3. While it is tempting to ascribe the lack of IFN production to the increased 
ISGylation seen in USP18"" cells, this is unlikely. Firstly, USP18 was shown to regulate type I 
IFN induction independent of isopeptidase activity and that ISG15 antiviral activity is independ¬ 
ent of IFN signalling (22, 43). Secondly, the observed lack of a type I IFN response may also be 
due to differences in the immune cell phenotype of USP18" " mice. Cong et al described a 
USP 18-dependent but ISGylation-independent difference in dendritic cell maturation (54). Simi¬ 
larly, we have observed a fundamental difference in the phenotype of USP18 _/ " macrophages and 
kupffer cells, with a shift from the Ml pro-inflammatory phenotype in USP18 +/+ mice to an M2 
regulatory phenotype in USP 18" mice (manuscript in preparation). This result may suggest that 
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the fact that MHV-3 does not induce IFN is not specific to the virus but the phenotype of the 
USP18 _/ " macrophages and dendritic cells, both of which are key producers of hepatic IFN (23, 
53, 55). Differences in immune cell phenotype and the consequent lack of IFN production in 
USP18 _/ " mice raise the possibility of confounding immune effects in our in vivo data. However, 
the lack of a type I IFN response in the USP18" " mice would be expected to promote viral repli¬ 
cation, but we observed the reverse effect. The lack of IFN signalling thus allows us to attribute 
our results to ISGylation as opposed to IFN signalling in USP18" " mice. 

The mechanism by which ISGylation exerts an antiviral activity against MHV-3 is un¬ 
clear. ISGylation can modulate viral production in one of two ways: by being conjugated to viral 
proteins in a manner that disrupts the viral life cycle, or by conjugating to host proteins that play 
a role in the viral life cycle or in the host innate immune response (11). In the context of the IS¬ 
Gylation of viral proteins, the conjugation of ISG15 to the Influenza A NS1 viral protein is anti¬ 
viral to influenza A infection in humans (39). ISGylation ofNSl prevents its association with 
importin-a and its subsequent translocation to the nucleus, ultimately inhibiting viral replication 
(38). In an example of host protein modification, the ISGylation of host 4E homologous protein 
(4EHP) protein exerts an antiviral effect by interfering with viral (and host) protein translation. 
4EHP is a 5' mRNA cap binding protein that competes with eukaryotic initiation factors-4 (eIF4) 
to suppress translation. The ISGylation of 4EHP increases its RNA binding capacity, and 
strongly suppresses translation of both host and viral proteins during the innate immune response 
(37). Lastly, another group reported a mechanism that incorporates both host and viral protein 
ISGylation: the antiviral pressure of ISGylation observed in this model is exerted by the indis¬ 
criminate ISGylation of newly synthesized proteins that follows IFN stimulation or USP18 dele¬ 
tion (36). The nonspecific ISGylation of the nascent subset of each protein population present in 
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the cell has been suggested to disrupt viral proliferation through a dominant negative effect. Such 
an effect has been shown in models of human papilloma vims (HPV) and HIV infection (36, 56) 
and is consistent with our finding that the effect of ISGylation blunts viral replication but does 
not abrogate it altogether. Further studies are necessary to determine the mechanism by which 
ISGylation specifically delays MHV-3 replication. 

Although ISGylation exerts an antiviral pressure on the MHV-3 coronavirus, the virus 
ultimately overcomes this pressure both in vitro and in vivo. Viral titers in states of increased 
ISGylation (USP18 _/ ~), approach peak levels seen in states of low ISGylation (USP18 +/+ ) follow¬ 
ing a 3-day delay. Additionally, N protein production and PLP2 expression levels increase while 
ISGylation levels in USP18 _/ " mice decrease (Figure 6-7). All Coronaviruses express a DUB that 
removes both ubiquitin and ISG15 from conjugated proteins (29, 57, 58). PLPs are multifunc¬ 
tional viral proteins and have several functions relevant to coronavirus production within the host 
cell including the processing of the viral polyprotein. In our study, some of the inhibitory effect 
on MHV-3 production that follows treatment with PLpro inhibitors may be due to this effect (34, 
51). However, treatment of wildtype PEM and hepatocytes with PLpro inhibitors increased cellu¬ 
lar ISGylation following infection, and inhibited viral N protein expression (Fig 5). Our data 
therefore, may be the result of both nsp3 cleavage and ISGylation but, it definitely contributes to 
the stripping of ISG15 from target proteins. Regardless of the relative contributions of each 
mode of action, these data support the use of PLP2 inhibitors as one arm of treatment against 
coronavirus infection. The DUB activity may thus contribute to the evasion of the host antiviral 
response by altering immune responses through delSGylation or by reactivating ISGylated viral 
proteins. 
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We have previously shown that increasing ubiquitination exerts an antiviral effect on 
coronaviras replication (28), and in this study we have shown that ISGylation exerts a similar 
effect. Both DUB functions of PLP2 would favour MHV3 replication, though the PLP2 has more 
affinity for ISGylated proteins than for ubiquitinated proteins and therefore may be more rele¬ 
vant as an antiviral mechanism (29, 30, 59). 

Overall our data are consistent with the ability of ISGylation to exert an antiviral pressure 
in a model of severe MHV3 coronaviras infection. This effect is most marked in the setting of 
USP18 deletion, suggesting that targeting USP18 might have a clinical benefit in treating severe 
Coronaviras infections. However, our data also strongly argue that treatment of a Coronaviras 
infection should be multi-faceted, since the virus has evolved mechanisms to overwhelm the an¬ 
tiviral effect of ISGylation - one of the chief effector mechanisms of IFN. As illustrated by this 
study, coronaviras can escape the host antiviral response in a number of ways that are independ¬ 
ent of IFN signalling. IFN therapy, the current treatment for Coronaviras infection, should be 
considered in conjunction with other therapies, such as the inhibition of PLP activity and protea- 
some inhibition to comprehensively counter the effects of the viral DUB (28). 
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Figure Legends 

Figure 1: USP18"" mice showed delayed morbidity and death following MHV-3 infection 

A) Survival curve of USP18 +/+ and USP18 _/ " mice following infection with 50 pfu. MHV-3 i.p.. 
USP 18"" survival was prolonged up to 10 days with a median survival of 8 days compared to 4 
days for control mice. n=10 per group. p<0.0001 using the Logrank test for trend. B) H&E stain¬ 
ing of livers from USP18 +/+ and USP18" /_ mice at days 0, 3, 4, and 7 p.i.. Black arrows point to 
areas of focal necrosis. Photomicrographs at day 4 and day 7 for USP + and USP"’ mice respec¬ 
tively show confluent necrosis of the liver. Magnification 40x. Photomicrographs are representa¬ 
tive of each study group. C&D) The increase in liver markers of injury, AST and ALT, was de¬ 
layed in USP 18" " mice. AST and ALT were measured and represented as log(U/L) + SD. E) 
USP18 _/ " mice show decreased viral titers. Viral titers were measured daily until day 4 and on 
day 7 only in USP18 _/ " mice using a standard plaque assay. Data is represented as the mean 
log(PFU/g) ± SD. *** p<0.001, ** p<0.01, *p<0.05. 

Figure 2: The IFN response to MFEV-3 is blunted in USP18" " mice 

Animals infected with MHV-3 were sacrificed on Day 0 through 4 with an additional measure¬ 
ment at day 7 for USP18 _/ " mice (USP18 +/+ succumbed to disease by day 4 p.i.). 

(A) mRNA expression of both type I and type II IFN is inhibited in USP18"" mice. mRNA was 
harvested from livers at USP18 +/+ or USP18" /_ mice infected with 50 PFU. of MHV-3 and sacri- 
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ficed day 3 p.i.. (B-D) Serum levels of IFN were measured by ELISA. (B) IFN-a, (C) IFN-(3 and 
(D) IFN-y levels in the plasma were absent in the USP18’ /_ mice whereas USP18 +/+ were able to 
induce type I and type II IFN responses. Data is represented as mean ± SD. (A) n=2 (B) n>2, (C) 
n>3, (D) n>3. *** p<0.001, ** p<0.01, *p<0.05. 

Figure 3: Elevated ISGylation in USP18coincides with increased resistance to MHV-3 
infection in vitro. 

(A) USP18’ " primary hepatocytes are more resistant to MHV-3 infection than USP18 +,+ mice. 
Primary hepatocytes were infected with MHV-3 (MOI=l) and viral titers were assayed using a 
standard plaque assay. Titering was represented by mean ± SD. n=3. (B) Elevated ISGylation 
coincides with decreased N protein production in USP18" " hepatocytes. Time course of N protein 
production and ISG15 conjugate formation by WB. Hepatocytes were inoculated with MOI=l 
and harvested at the indicated time points. N protein production is a marker of viral replication. 
WB is a representative blot of two independent experiments. *** p<0.001, ** p<0.01, *p<0.05. 

Figure 4: ISGylation inhibits MHV-3 replication 

The El conjugation enzyme, UbelL, is critical for ISG15 conjugation to target proteins. UbelL 
expression was silenced to determine if increased expression of ISG15 or increased ISGylation is 
antiviral to coronavims. Primary hepatocytes were transfected with UbelL or control siRNA 
(Santa Cruz) and incubated for 24h. Hepatocytes were infected with MHV-3 (MOI=l) and cells 
were harvested 12h p.i.. Samples were visualized on a WB for UbelL expression, N protein pro¬ 
duction and ISGylation. Blot is representative of two independent experiments. 
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Figure 5: PLP2 inhibitors and ISG15 antiviral activity inhibit MHV-3 N protein produc¬ 
tion. 

PEM from USP18 +/+ or USP18" 7 " mice were inoculated with MHV-3 (MOI=0.1) for lh, the virus 
was then washed away and PEM were treated with PLP2 inhibitors (1 OOuM inhibitor 6, 7 or 
GRL0617) or IFN-a (lOOIU/ml) (34). Cells were harvested 9h p.i.. N protein and ISG15 conju¬ 
gates were visualized by WB. 

Figure 6: ISGylation delays MHV-3 replication 

Increased ISGylation coincides with prolonged survival and delayed N protein accumulation. 
USP18 +/+ or USP18" 7 ’ mice were infected with 50 PFU MHV-3/mouse. Mice (n=4) from each 
group were sacrificed daily and livers were harvested for WB analysis. N protein and ISG15 con¬ 
jugates were visualized by WB. 

Figure 7: The onset of PLP2 expression is delayed in USP18' ' mice 

N protein expression coincides with increased viral expression. To determine if PLP2 expression 
followed this expression pattern we tested mRNA expression of PLP2 following MHV-3 expres¬ 
sion. USP18 +/+ or USP18" 7 " mice (n=4/group) were infected with MHV-3 (50 PFU MHV- 
3/mouse). Livers were harvested on the indicted days and mRNA was extracted. Data is repre¬ 
sented as the fold increase of PLP2 mRNA over the housekeeping gene GAPDH. *** p<0.001. 
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